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Abstract 

We consider the unitary matrix model in the limit where the size of the matrices become infinite and 
in the critical situation when a new spectral band is about to emerge. In previous works the number 
of expected eigenvalues in a neighborhood of the band was fixed and finite, a situation that was termed 
"birth of a cut" or "first colonization" . We now consider the transitional regime where this microscopic 
population in the new band grows without bounds but at a slower rate than the size of the matrix. 
The local population in the new band organizes in a "mesoscopic" regime, in between the macroscopic 
behavior of the full system and the previously studied microscopic one. The mesoscopic colony may 
form a finite number of new bands, with a maximum number dictated by the degree of criticality of the 
original potential. We describe the delicate scaling limit that realizes/controls the mesoscopic colony. 
The method we use is the steepest descent analysis of the Riemann-Hilbert problem that is satisfied by 
the associated orthogonal polynomials. 
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4 Outer and local parametrices near the turning points 

4.1 Universality of the behavior 



5 Conclusion and generalizations [14] 
A Differential equations for the endpoints [lfs] 

1 Introduction and result 

The phenomenon that we want to investigate in this paper goes under the name of "birth of a cut" 
[101 1161 [5] or "colonization of an outpost" [5J [3] , namely the transition when one or more new spectral 
bands open in the asymptotic spectrum of the model. In particular wc want to focus on the transition 
between the microscopic regime (of finite number of eigenvalues) and the macroscopic regime (where the 
number of eigenvalues scales like N); we call the intermediate regime the mesoscopic regime. 

While the paper does not aim at being propaedeutic to the topic of random matrices, in this section 
we recall some general facts about the unitary random matrix model so as to set the context. Unitary 
random matrix model is defined by the probability distribution 



Z-\f cxp (-^TiV(M)j dM, Z n>N = J cxp (-^TrV{M)j dM , 



(1.1) 



on the space Tt n of Hermitian n x n matrices M, with V a real analytic function (the potential) that 
satisfies 

V(x) 

hm - — — — = +00. 

x^±oo [Qg[x z + 1) 

The eigenvalues x\, . . . ,x n of the matrices in this ensemble are distributed according to the probability 
distribution (See, e.g. [15], [9]) 

V^ N Hx u x n )d n x = Z-y-^U ^(x«) Y[( Xj - x k ) 2 d Xl . . . dx n , (1.2) 

j<k 

where Z„jv is the normalization constant. 

The correlation functions of the eigenvalues are related to orthogonal polynomials (see e.g. [9], [15]): 
let {7r„(x)}„ e N be the degree n monic orthogonal polynomials with weight e"" l/ ' 1 ' on R. [19] 

■Kn{x)-n: m {x)er ~ v ^ dx = h n 5 nm . (1.3) 

Let us construct the correlation kernel by 

K ntN {x,x') = e -W<yW +v V» J2 nj{x) * j{x '\ 

3=0 3 
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By the Christoffel-Darboux formula, this kernel can be expressed in terms of the two orthogonal polyno- 
mials 7r„(x) and Tr n -i(x) instead of the whole sum: 



N 



r x x '\ = c -2 ! ^(V{x) + V(x')) 7T n{x)TTn-l(x') - TT n {x')'K n -l(x) 

h n _i(x - x') 



(1.4) 



The basis of our analysis relies on the Fokas-Its-Kitaev formulation [ITJ [12] of OPs in terms of the 
following RHP for the 2x2 matrix Y(z) (for brevity, we drop the explicit dependence of Y on n) 



y+(x)=Y_0r) 



1 e ~§V(x) 

1 



z n 
z~ n 



(1.5) 



(1.6) 



and the polynomial 7T n (z) is simply Yn(z), while the kernel is recovered from 

— ZTT1(X — X ) 

Then the m-point joint probability distribution function can be written as the determinant of the kernel 

T^m N \ X ^- • • ,^m) ■= det (K^Nix^Xk))^^^ 



In the limit lim„ 



n,N — >oo 



77 = 1, the eigenvalue density 



■R.\ n -">(x 



of the ensemble (|1.1[) is asymptotic to 



the equilibrium measure p(x) [8], |13j . |18j 



Um 22 )El = p (a.) 

.,JV->oo,£— ►! n 



where the p{x)dx = dp m i n {x) is the normalized density of the unique measure p m in(x) that minimizes 
the energy 

I(/i)=-T / / log |z- 2/1^(30^(2/)+ / ^(a;)^^) 



among all Borel probability measures p on R. The fact that p m in{x) admits a probability density follows 
from the assumption that V{x) is real and analytic [6]. Moreover, it was shown ibidem that for real and 
analytic V(x), the equilibrium measure is supported on a finite union of intervals. 

1.1 Colonization at an outpost 

The following conditions are satisfied by the equilibrium density p(x) [5] , [18] 

2T / log \x - s\p(s)ds - V(x) =£, x e Supp(p), 

7 (1.7) 

2T / log \x - s\p(s)ds - V(x) <£, x e M/Supp(p), 

Jg. 
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for some constant I (also known as Robin's constant). For a generic potential V(x), the inequality in 
(|1.7p is satisfied strictly outside the support. Suppose however that there is some point xq S\vpp(p(x)) 
where the inequality is not strict 

2T / log \x - s\p(s)ds - V(x ) = I. (1.8) 

Such a potential V is called irregular ([7]); a small perturbation of the potential may induce a new 
interval of support of p to form around xo- We may think of this phenomenon as the eigenvalues 
colonizing the point xq, which we will call the outpost. This situation has been considered previously 
and the term 'birth of new cut' was used in some of the studies [5], [TU], [To] . 

In the studies j0, PH], [IS], the colonization phenomenon was considered when a finite number 

of eigenvalues start appearing in the outpost xq- It was shown that the eigenvalue statistics near the 
outpost can be described by that of a finite size Hermitian matrix ensemble, or a microscopic ensemble. 

1.2 Genus transition in random matrix models: the proliferation of a colony 

Some interesting questions about how new intervals in the 
support are forming remained unanswered. For example, 
whether several intervals in the support can form simultane- 
ously or they have to form one after another; how to describe 
the eigenvalue statistics when the number of eigenvalues in 
the colony becomes large. This paper aims at addressing 
some of these questions, namely, we want to analyze the tran- 
sition from the 'first colonization' to the situation where one 
or more new intervals are fully formed. A schematic view is 
shown in Fig. [TJ Since the size of the colony, though small 
compared to the main cut, is taken o(N) and unbounded, we 
use the term mesoscopic colonization. 

The setup is as follows: let V{x) be a critical potential 
such that (|1.8| is satisfied at a point xq outside of the support 

of the equilibrium measure and let the order of vanishing of 

iigure 1: Mesoscopic colonization. 

(|1.8|) be 2v + 2. In particular, at the outpost, we have 

<j>(x) = V(x) -2T J log \x - s\p(s)ds + I = C Q (x - x f u+2 {l + 0((x - x )) , C > . (1.9) 

Without loss of generality we will perform a translation of the problem so that xq = 0, but we will keep 
referring to it at xa not to confuse it with other zeroes. The function <j)(x) is called effective potential 
since it represents the sum of the external potential V and the Coulomb (two-dimensional) potential 
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generated by the equilibrium distribution. Let Bj(i) be a bump function that is 1 inside an interval 
J C K around xq and outside an interval J D J around xq. Both J, J are chosen small enough so as 
not to contain any point of the support of p. We will study perturbed model 



(1.10) 



Z~n exp [~Tr&(M)J dM, Z n , N = cxp [~ TrV(M)) dM, 

where V(x) is a 1-parameter perturbation of V(x) (See Figured]) 

V(x) =V(x)+Bj(x)A K>N (x), (1.11) 

where k is of order 0(N 1 ^ t ) with < t < 1, and A(x) is analytic near x = xq. Due to (| 1 .9[) . we can 
define a local parameter r\ inside a finite neighborhood D around xq as follows 

V= (S 0(X) ) 7 = ^ + (2:eID)) ' (L12) 

kT \ 7 1 . , 

7:= 7^^- (1.13) 



C N J ' 2i/ 

We will show that for a suitable choice of the perturbation function A{x) the eigenvalues of the matrix 
model (| 1 . 1 0(1 arc distributed on micro-cuts in D whose image on the 77-plane is a collection of at most 
v + 1 segments. In terms of the coordinate x, this support shrinks at a rate 0((k/-/V) 7 ) = 0(N~ tl ). 

The bump function Bj is used to keep the technicalities to its minimum and is not essential to the 
construction: changing between two such bump functions will introduce a difference in the description 
which is exponentially small (as N — > 00) and hence beyond all orders of perturbation. Such a manipu- 
lation of the potential is very useful in handling the otherwise complicated "double scaling limit" . The 
result of this paper can be encompassed in the following theorem: 

Theorem 1.1. Let V be real analytic and irregular, with an effective potential vanishing at xq as dictated 
in \1.9\) and let k = kn be a sequence of integer^ such that n = 0(N 1 ~ t ) for some 1 > t > 0. Let Tj be 
the scaling coordinate near xq given by eq. WTty and V mes (rj) = n 2 "+ 2 + J2%Y hrf > h e R be a real 
monic polynomial potential of degree 2v + 2. Let fi m es{v) be its equilibrium measure minimizing 

V m es(v)dp m es(v) ~ log |77 - ■q l \dHmes{'n)dHmes{'n l ) ■ C 1 - 14 ) 

J JRxR 

Let D be a small neighborhood of Xq in the complex plane on which i] (|1.12[) is conformal. Then we can 
define a function A(x) = Ajv l/t (x) as in Def. \3.2\ to be used in (1.11)) that is analytic and bounded on 
J, of order 0(e) = 0((k/N) 7 ) = 0(N^ 7t ) (and uniformly so w.r.t. N,k and x £ J) such that in the 
limit N — > 00, k = 0(N 1 ~ t ), < t < 1, the eigenvalue density of the matrix model (1.10)) in D is 
asymptotically given by the equilibrium measure fJ-mesiv) * n ^ e n-plane. 



3 The requirement is purely technical and could be disposed of, at the price of complicating the analysis, without changing 
the result. 
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The proof takes up the whole paper. The explicit form of A(x) is given later in (|3.14D . 

We will argue in Sec. 14.11 that the -admittedly contrived- construction of the perturbation is in fact 
sufficient to capture the universal behavior. 

The eigenvalues of (|1.10p . except the ones on the macroscopic cuts, are on the support of dp mes {rf). 
From (|1.12p . we see that the lengths of these cuts are of order (^r) 7 - 

Remark 1.1. Theorem \1.1\ means that for a suitably chosen perturbation of the critical potential, the 
number of micro-cuts that are formed depends on the mesoscopic potential V mes (x) . In particular, multiple 
cuts can be formed simultaneously if the equilibrium measure of V mes (x) is supported on multiple cuts. 

2 Equilibrium measure in the mesoscopic problem 

Consider the mesoscopic potential in Theorem 11.11 with F(rf) := X^j^i tjrf > 

V mcs (ij) = n 2 »+ 2 + F(n) , degFM <2v+l . (2.1) 
We define as usual the corresponding g-function as the logarithmic transform of the equilibrium measure: 

ffmes(C) = / ln(£ - ?7)Pmcs(??)d?7 (2.2) 



where p m es(j]) is the probability measure on K that minimizes the familiar energy functional 

Fmes[p\~ J V mcs (rj)p mcs (7])d7] + J J m T~ g| Pmcs(?7)|Pmcs(£) d?7 d£ (2.3) 

The support of p mcs is a finite union of intervals [5] and it is possible to see that in fact it can have 
at most v + 1 disjoint intervals. The g-function has an expansion for large argument of the form 

9mes(0 = k£ " E H ' h i := / P™Mr? dv- (2.4) 

j=X 

The mesoscopic equilibrium measure satisfies the same inequalities as in (| 1 . 7|) with an appropriate (meso- 
scopic Robin's) constant £ me s' 

V mes (r)) - 2$tg mcs (n) + £ mcs = , 77 <E Supp(p mcs ) , (2.5) 
V mes (r)) - 2Kg mcs (77) + £ mC s > , ?y ^ Supp(p mcs ) . (2.6) 

We will need the following truncation of the expansion 

Definition 2.1. The truncated mesoscopic g-function is defined as 



J 



Note that we have defined both g mes and the function f m es ■ 

The minimal level of truncation k will be determined in (|3.23p . but for the time being it is a parameter 
of our problem. 
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3 Singularly perturbed variational problem 



In order to construct the deformation of the original problem so that we obtain the desired double-scaling 
limit we need to work a bit more compared to [2j [3] . In particular the global (/-function will be modified 
to a certain extent because mesoscopic colony is "big" enough to affect the minimization problem for the 
macroscopic spectrum. 

Let V(x) be a real-analytic potential. It is known from [5] that the support of the corresponding 
equilibrium measure consists of a finite union of disjoint finite intervals |J[q;2j-i , 012 j]- We define the 
complexified effective potential by the formula 

ip{z) := V{z) -2T j p(t) ln(z -t)dt + £ (3.1) 
Jr 

Due to the multivalucdness of the logarithm tp is only defined on a simply connected domain, which 
customarily is chosen as C \ (—00, maxsupp (p)] [7J. If the point xq belongs to a finite spectral gap (a, b) 
(a finite connected component of M \ supp (p)) we can alternatively define <p as a holomorphic function 
on C \ {(— 00, a] U [6,00)}; the only effect in this re-definition is to modify the so-called model problem 
(or outer parametrix) by a constant (in z) multiplier. It is more convenient for the discussion (but not 
at all crucial) to use a definition of <p which is analytic at xq and so we will assume this is the case. For 
example, if xq lies to the right of supp p then no additional complication arises. 

The condition that V(x) is irregular (|1.9p at xq $ supp (p) is translated in terms of ip as the condition 
p(x) = C x 2v+2 (l + 0(x)), Co > 0. 

We will assume, for simplicity, that all other turning points are simple, namely at the endpoints of 
the intervals of the support of p we have <p'(x) ~ C(x — oij)^ (1 + 0(x — ay)). 

The goal of this section is to define a small perturbation to the unperturbed (/-function (or the 
unperturbed effective potential) which will serve to normalize eventually- the RHP for the corresponding 
orthogonal polynomials. 



Definition 3.1. The mesoscopic conformal scaling parameter rj is defined by the following equations 

IV / N \ 7 r 

W +2 = 7^0*0 V ■= = ^(x) = -(1 + 0(x)), (3.2) 



7 := 2^ ' 6 := (£n) • (3 - 3) 

The choice of symbol (rj) is made on purpose to match the use of coordinate that was made in the 
previous section (Sec. [2]). Define the following Laurent polynomial in x 

f{x/e) := _ res /mcs(??(z)) dz = J2 T% . fnMv) as in m (3.4) 
z=o z — x * — ' (x e)J 

j=i 

We note that f3j = bj + O(e) are analytic functions near e = 0. 
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The singularly perturbed minimization problem consists now in minimizing the following functional 

=:V(x) 



T t := J (v(t)-2^H e (x)) d^(t)+T J J d»(t) dfx(s) In ^ 
H e (x) :=ln\x/e\-f(x/e) , 

d/i(t) = 1 - , supp( y u) Cl\ J. (3.5) 

where the minimization is taken over the set of Borel measure that is supported on R \ J. Note that, 
with the above definition of H e (x), the following property is verified. 

H e (x) - g mcs (v) = O(e) + 0(x)(l + 0(e)) . (3.6) 

We point out that the potentials V = V(x, e) are admissible on R \ J in the sense of potential theory 
f|18|) for sufficiently small e: let p e be the corresponding equilibrium measures. Then we can define the 
modified g-function by 

g{x) := [ \og{x -t) Pe (t)dt , (3.7) 

JR 

and the modified effective potential by 

<p(x) := V{x) - 2Tg(x) + £ + T^£ mcs (3.8) 

where we have written the Robin constant for the modified minimization problem as £ + T-^^ mcs for 
convenience, and -by definition- it is such that its real part of tp is zero on supp(p e )- Note that for 
K = (e = 0) , the solution of the variational problem (|3.5[) and the original one over the whole real axis 
coincide since both fulfill cqs. (|1.7p . 

We can then apply the results of [14] to conclude that V is a regular potential (for e small) on R \ J. In 
particular we quote the relevant 

Theorem 3.1 ( [14] . Theorem 1.2). Suppose V and V n , n = 1,2,... are real analytic external fields on 
R such that the following hold: 

1. V n and the first three derivatives converge to V uniformly on compact subsets o/R; 

2. The growth condition limi^i^oo V n (x)/ In \x\ = +oo holds uniformly in n. 

Then the supports of the corresponding equilibrium measures are uniformly bounded. Furthermore, if V 
is regular then so are all the V n eventually. 

There are two remarks due at this point 
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• Clearly we can adapt the above theorem to any family of potentials V e with the obvious modifications 
of the statement; in this case an isotopy argument implies that if Vb is regular and has K component 
in the support of the corresponding equilibrium measure, so happens for V e , as long as e is sufficiently 
small. 

• Thm. 13.11 is stated on R but, reading the proof in [14| . it appears that there is no difficulty in 
replacing R with R \ J. or even any union of intervals, for what matters. 

Thm. 13.11 with the above trivial extensions implies that the number of components of the supports for p e 
is finite and constant (for e sufficiently small) and the endpoints are smooth functions of e. In fact it is 
possible to derive (nonlinear) differential equations for the endpoints as functions of e. In the appendix 
we give the result without proof, since it is not necessary to the considerations to follow. 

3.1 Modified orthogonal polynomials 

We choose a small interval J around the outpost that does not contain any other cndpoint. We will 
consider the following modified orthogonality relations 

h„S nm = / p n (x)p m {x)e-T rV{x) dx (3.9) 

where the perturbed potential V(x) was given in (|1.11|) . 

For simplicity we will also assume that k = kjv depends on N in such a way that 

• K = /tjv is an integer; 

• k n = OiN 1 -*), 1 > t > 0. 

Were we to allow k to be non-integer, we would have to complicate the analysis by taking into account 
that when k crosses the half-integers an improved local parametrix needs to be used as in [5] . This would 
only lengthen (considerably) the paper while providing no further insight into the phenomenon we want 
to describe. 

3.2 Dressing the RHP with the singularly perturbed g— function 

For the orthogonal polynomials at ()3.9|) we take the RHP (|1.5p for Y with V instead of V. 
We define 

:= e-^e- 24 *" ^-" * Y{z) ^{b^)-^)^-^,)-^ (3 1Q) 

The various prefactors of Y(z) above are only to ensure that &(z) = 1 + 0{z~ l ) 1 the e~ K<T3 term coming 
to compensate the term ln(x/e) that appears in H e (j3.5[) . In this way, the ^-functions is "stripping off" 
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the outer paramctrix from "all" the zeros including the ones at the outpost. This approach is different 
from the one in [3J [3] and actually closer to [HUE]. As a result the jumps on J for ^(z) become 



*+(x) = *_(x) 



1 e -T {V-2Tg+I)+K(2H e (x)-e meB ) 

1 



(3.11) 



By virtue of the variational problem that g solves and since V = V outside of J, the analysis on the 
support of p can be carried out in verbatim as in [7] , keeping in mind that the cndpoints are slowly varying 
functions of k/N. We want to focus on the problem near the outpost, as it contains the whole essence of 
the new phenomenon. 

In order to have locally the (simplest form of the) RHP for the mcsoscopic potential we need to have 



(V(x) - 2Tg(x) +Tj- 2kH € {x) + K £ mes = k (V mes {r)) - 2^0?) + (3.12) 



N 
f 

N 



— {V(x) - 2Tg{x) + t) + K F{ V ) - 2Kg mes (v) + Ki mes . 

V v - 

— KT} 2lJ + 2 

(3.13) 

Simplifying the above expression, we have 

Definition 3.2. The perturbation function A(z) := A k .n(z) to be used in the perturbed potential 111.11]) 
is defined as 

A(z) = V — V = — (2T(g - g) + (£-£))+ ^ (f( v ) + 2 (H e (x) - g m es(v)) ) (3.14) 

Remark 3.1. We recall what are the input data in Def. \3.2l so as to make clear the definition is not 
circular: we need 

• the unperturbed nonregular potential V(z) (with the property jl.9]) ): 

• the mesoscopic potential V mes (n) = rj 2v+2 + F(i]), with F(n) an arbitrarily chosen polynomial of 
2u+l; 



• the order of truncation k. 

The other functions appearing in {3.1$ are H e (defined in i3. 5]) ). the truncated g-function (Def. \2.1\) 

In order to manifest the analytic properties of A{x) = V — V we point out that the singularities of 
H e (z) and p mes cancel out precisely by (|3.6|) to give a locally analytic function in the neighborhood of 
the outpost. Also the largest deviation is given by the term ^F(ri) = e 2lJ+2 0{e~ 2l '~ 1 ) = 0(e). The 
deviation 2T(g — g) + (£ — £) is bounded by 0(e 2v+2 ) in general (see, for instance, the lemma lA - !]) . 



10 



With this position for A we have the following RHP in the neighborhood of the outpost 

I e -K(Vmes(?;)-2gmes(r;)+^mes) 
1 

tf (z) = O{z )c- KH ' (z)a3 = 0{if )c- K ° mcs{ri),J3 , 77 -> . (3.15) 



The growth behavior at the origin is obtained from the definition of 'J (|3.10p . The reason why the two 
behaviors at z = on the second line of (|3 . 1 5[) are equivalent is due to the fact that H t (z)\s precisely the 
singular part of g mes , as follows from (|3. 5[) and (|3.4p . 



3.3 Local parametrix at the mesoscopic colony 

Let Pj(rj) = Pj(fj', k) be the monic orthogonal polynomials for the (varying) measure e _KVmos( '' , ^ drj. We 
want to construct an exact solution R K (r]) to the jump condition (|3 . 1 5|) such that on <9B it behaves as 
1 + 0{N~ a ) for some positive a and uniformly in 77. 
Consider the matrix 



R K (r,) = e"^» 



P K ( V ) C[P K ]{rj) 



h K -i K ~ LVIJ h K -i 

It is immediate to verify that it solves the following jump condition and asymptotic behavior 

■y c -K(Vmes(r?)-2gmes(rf)+4i 



R K { V )+ - R K ( V )- 



1 



(3.17) 



R K {r]) = l + o(^\ 7] ^00. (3.18) 

R K ( V ) = 0(l)e- KS ™*° 3 , 77 -0, (3.19) 

The last error term C(7/ _1 ) in (|3.18p depends on k and we need to control this uniformly in the limit 
K — > 00. To do this, which is a crucial fact to the error analysis of the asymptotics, we use the standard 
knowledge on the asymptotic behavior of the mesoscopic orthogonal polynomials, that comes from the 
steepest descent analysis [7] of the local RHP with respect to the local coordinate 77. Briefly this amount 
to say that, as k — > 00 wc can obtain expressions 

R K (ri) =(l + f e^e^mcsM^e-KSmcsM^ ?? ^ (3 20) 

Note that mesoscopic Robin's constant disappears by virtue of our well-crafted choice of perturbation. 

The first factor in (|3.20p comes from the error of the mesoscopic error matrix and E mos is determined 
by the nature of the mesoscopic system; for a usual situation with regular mesoscopic potential, we have 
E m cs = 1- The 2x2 matrix here denoted by 6(77) is the theta function expression for the asymptotic 
of the orthogonal polynomials that solves the so-called "model problem" , with jumps on the support of 
fi m es and in the interval between (See [7] or [3]). It is known that the factor 8(77) behaves as 

9(77) = 1 + e>(77- 1 ), 77^00. 
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where the error term is also bounded in k as k — > oo. 

The trailing exponential factors in (|3.20[) determine the minimal order of the truncation (|2.7p : using 
Definition 12. 1[ we have g mea (rf) — g m cs(v) = 0(l/r/ k+1 ). Therefore, for r\ on the boundary SB we have 
the uniform estimate 

R K (r 1 )=(l + o(-^-)) (l + o(-\ 1(1 + 



T) ) ) \ \7] k+1 



= l+O — - l + O 



( 

7 (fc+l)+l 

l + O' 



7V7(fe+i) 



(3.21) 



V i / 

The last contribution to the error term marked with ★ in (|3.21[) is the most important one: demanding 
that the error decays imposes a condition on the minimal k of the truncation in Def. (|2.ip . Indeed, the 
growth of k must be such that there exists a minimal fc m m for which the last term is o(l). In other words 
the order of growth of the colony must be 

K < iV 1_ Wd+i . (3.22) 

for k sufficiently large. 

If k = 0(N 1 ~ t ) for some < t < 1 then we need to choose k so that 

^>+Tyn <f ° ^(^(H -1 (3 - 23) 

This determines the minimal order of truncation in (|2.1[) and in all the analysis that followed. The 
error bound is then dominated by the last term in (|3.2ip and can be made as close as desired to 0(k 7 /N 1 ) 
by choosing k sufficiently large. For instance, choosing the next-to-minimal k -which we do henceforth 
yields an error bound 1 + 0(iV~ 7 *). 

4 Outer and local parametrices near the turning points 

We will not go into much detail regarding the rest of the asymptotic analysis because it is quite standard. 

In fact the strong asymptotic for ^ in (|3.10p is obtained in the identical way as in [7] by "opening 
the lenses" around the intervals of the support for the perturbed variational problem (|3.5p . As we have 
assumed, the variational problem is regular outside of a the (^-neighborhood of the outpost, hence the 
procedure is verbatim as in [7] . The only caveat is that the endpoints are slowly varying functions of the 
small parameter = 0{N~ r ). 

It should be clear that the error term of the analysis becomes 0(1) as t — > 0+, namely, as k grows 
at the same order as N (at which point the new gaps must be "fully formed"): more and more terms 
need to be added to the truncation 12. H Eventually one must solve an exact minimization problem 
when the colony is fully grown and the transition from mesoscopic to macroscopic will be complete. 
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The outer parametrix i , / *j—2*//(<y(2k+l)+i)\ 

is the same as for s~~~v' V ^ 



OiN- 1 ) 

j r\ 

the unperturbed prob- y \ -X—-' ^ \ / 

leni [7], as well as all \ X; \ / 

the local parametrices ~ \\ / 

near the endpoints of l + 0(iV _00 ) 

the spectrum. 

In particular. de- 



Figure 2: The order of the jumps of typical residual Riemann Hilbcrt problem 
for the error term, here depicted for a one-cut situation with regular endpoints 
noting by ^ the outer (where the Airy local parametrix can be used). 

parametrix, the global 
parametrix will be of the form 

!^>°°(z) outside the local regions 
^°°(z)A(z) near the endpoints , (4.1) 

^°°(z)R K (z) at the outpost 

where A(x) is the parametrix constructed by Airy functions as in [7j. Here is the outer global 
parametrix constructed out of theta functions as in [7] or equivalently using spinors as in [1] . 

The error matrix £ := \f» as vl/ -1 has several residual jumps: in particular on the boundary of the disk 
at the outpost we have 

S+iS-Y 1 = * as i? K (* as ) _1 = 1 + 0(N-^) . (4.2) 

We now comment on why the choice of bump function is totally irrelevant; indeed, on the real axis 
and outside of D the jump is exponentially close to the identity matrix (uniformly). Changing bump 
function trades such jump by another one, equally close to the identity jump, while leaving the jumps 
identical within B. As a consequence, the ratio of the solutions corresponding to different choices of 
bump functions would solve a RHP with jumps exponentially close to the identity everywhere (and in 
L°° n L 2 ). Thus the two solutions would differ only by exponentially suppressed terms, well beyond any 
order of perturbation. 

4.1 Universality of the behavior 

The perturbation of the potential has been chosen in the contrived form (jl.lll Dcf. 13. 2jl to eventually 
yield the simplest form for the local Riemann-Hilbert problem in Sec. 13.31 the gist of all the construction 
is such that in the scaling coordinate 77 the jump on the interval J is given precisely by (|3.17|) . As 
often happens [2l [3], the logic of our construction is slightly backwards from the more conventional 
approaches |10[ 1161 E] : we guess what local RHP would give the the phenomenon we expect on heuristic 
grounds, and try to "reverse-engineer" the appropriate deformation of the potential. This approach, while 
completely rigorous and also simpler to implement, is possibly not the most transparent to the reader. 
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The perturbation A(z) (|3.2j) is (a) analytic in z and (b) of order (k/N) 2 "+ 2 : from a heuristic point of 
view (based also on similar setups in the study of the universal unfolding of singularities [1]) it is natural 
to expect that these are the only relevant features to generate a mesoscopic colonization. Of course there 
is much more detailed information that goes into our approach, because the mesoscopic colonization 
as we described — with a fixed (i.e. non-scaling in N) local matrix model — can only be obtained as a 
multi-scaling limit; isolating k eigenvalues requires one scaling, forming a specific local cut structure will 
require rather complicated scalings. While we could not find a simpler, more direct path that starts 
from the perturbation and ends at a full description of the scaling regime, we do not expect that such a 
description, while logically more appealing, would be any simpler. 

5 Conclusion and generalizations 

• A quite parallel analysis could be performed in the case of the colonization of a hard-edge as in [3] . 
While the logic is identical, there are sufficient small details that would require a separate analysis, 
but with the final picture being completely analogous: in that case too one can have -depending 
on the degree of irregularity of the unperturbed potential- a mesoscopic growth of several meso- 
intervals for the equilibrium measure. We believe that the analysis is not sufficiently different to 
require a separate paper and yet not similar enough to put it here at the expense of clarity and 
conciseness. 

• It was also pointed out to one of u:jf| that the technique of analyzing the colonization (microscopic 
and mesoscopic) can be applied almost verbatim to the study of the trailing-cdgc of the solution of 
the Korteweg-deVries equation after the time of gradient catastrophe; we reserve to come back on 
this issue on a subsequent publication. 

• Since the mesoscopic potential can be an arbitrary polynomial, we could choose Vmes as a nonregular 
potential such that it has a point outside the support of pmes where the variational inequalities 
(|1.7p fail. Thus, one may have the whole picture of microscopic/mcsoscopic colonization within 
the analysis of the mesoscopic paramctrix. By perturbing Vrnes accordingly one could study a 
multiscale colonization. Since a polynomial potential of degree 2d+2 can have such a nonregular 
point with order at most 2d, we can "embed" the micro/mesoscopic pictures one into another at 
several nested scales at most 2v times, if the macroscopic potential has an irregular point as the 
one studied in this paper. We could call this multiscale situation the "MatryoshkeQ colonization" . 

4 We thank B. Dubrovin and T. Grava for the indication. 

5 A Matryoshka doll is Russian toy consisting in a set of dolls of decreasing sizes placed one inside the other. 
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A Differential equations for the endpoints 



Here we simply state a result (Prop. lA.ip that can be proved along the lines of the Buyarov-Rakhmanov 
equation We use the same notation as in the text (|3.8|) and we set 

y = 1 -n x) -TjeM^ = ^ x) . (L1) 

It is known that y solves a (pseudo) algebraic equation of the form 

2g+2 

y 2 = F e (x) 2 JJ(x- ttj -( e )) (1.2) 

3=1 

and F e (x) is a real-analytic function (depending on e) with a pole of degree k at x — 0. 

Lemma A.l. For small nonnegative e we have (the dot means differentiation in e) 

uj(x)dx:=ydx (1-3) 

is the unique meromorphic differential (whose existence follows from standard arguments in algebraic 
geometry) on the hyperelliptic Riemann surface branched at the endpoints w 2 = Ylj^i i x ~ a j) with the 
properties that 

1. all periods are purely imaginary; 

2. uj(x)dx has poles only at the point above x = with residues =f2T (respectively, on each sheet). 

3. at x = (on the physical sheet) it behaves as 

"(x) ~ ~ + -0^d t d x (f(x/c)) (1-4) 

Note that the second part contains poles of order strictly higher than one and hence corresponds to 
a second-kind differential. 



In particular we note that u>(x) can be written as 



l 3 

with R(x) a rational function of the form 



^ = ' 2„+2 = (L5) 



R(x) = Pfe+l ( - ) + P g -i(x) (1.6) 



and P m (Z) denotes some polynomial of degree m of the indeterminate Z. The above three facts completely 
determine R(x) as a function of 's, f3j 's and e. 
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Proposition A.l. The endpoints solve the following differential equation 



(1.7) 



F e (x) 



R(x) - F t (x) J2j Oj Uk^( x - 



(1.8) 
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